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ABSTRACT 

We present a series of RXTE observations of the nearby obscured Seyfert galaxies 
ESO103-G35, IC5063, NGC4507 and NGC7172. The period of monitoring ranges from 
seven days for NGC7172 up to about seven months for ESO103-G035. The spectra of 
all galaxies are well fit with a highly obscured (Njj > 10 23 cm -2 ) power-law and an 
Fe line at 6.4 keV. We find strong evidence for the presence of a reflection component 
in ESO103-G35 and NGC4507. The observed flux presents strong variability on day 
timescales in all objects. Spectral variability is also detected in the sense that the 
spectrum steepens with increasing flux similar to the behaviour witnessed in some 
Seyfert- 1 galaxies. 

Key words: . galaxies:active - X-ray:galaxies - galaxies: Seyfert - galaxies individual: 
ESO103-G35: IC5063: NGC4507; NGC7172 



1 INTRODUCTION 



Recent X-ray missions have brought a rapid progress in our 
understanding of the X-ray properties of Seyfert galaxies 
(for a review see Mushotzky, Done & Pounds 1993). The 
X-ray spectra of Seyfert 2 galaxies as observed by Ginga , 
ASCA and recently BeppoSAX have proved to be very com- 
plex (Smith & Done 1996; Turner et al. 1997). In broad 
terms most Seyfert 2 X-ray spectra can be well fitted by 
a power-law continuum (typically T ~ 1.8), plus an Fe-K 
emission line at 6.4 keV and a reflection component (e.g. 
Lightman & White 1988, George & Fabian 1991). This latter 
component, which suggests the presence of a large amount 
of cold material in the vicinity of the nucleus, flattens the 
observed continuum and can dominate the spectrum above 
~ 10 keV. In most Seyfert 2s the above emission compo- 
nents are viewed through a large absorbing column density, 
typically Nh > 10 23 cm -2 . This screen, which suppresses 
the soft X-ray emission through the process of photoelec- 
tric absorption in cool atomic and molecular gas, is possibly 
associated with a large (pc scale) molecular torus. In the 
limit of very high column densities, Thomson scattering will 
also diminish the more penetrating hard X-ray emission. In 
some sources additional emission in the form of a soft X- 
ray excess is observed below ~ 3 keV probably as a result of 
scattering of the intrinsic power-law continuum by a strongly 
photoionised medium (eg Griffiths et al. 1998 in the case of 
Markarian 3) . In order to observe such soft X-ray emission it 
is clearly a requirement that the scattering medium should 



extend beyond the bounds of the obscuration of the molecu- 
lar torus. The X-ray spectra of Seyfert-2 galaxies gave wide 
support to the standard AGN unification model (Antonucci 
& Miller 1985) . According to this paradigm, the nucleus (su- 
permassive black hole, accretion disc and broad line region) 
has basically the same structure in both type 1 and type 2 
Seyferts, but depending on the circumstances, can be hid- 
den from view by the molecular torus (Krolik & Begelman 
1986). Specifically if the source is observed at a sufficiently 
high inclination angle, and thus the line of sight intersects 
the torus, it would be classified as a Seyfert 2, whereas for 
all other orientations it would be deemed to be a Seyfert 1. 

In contrast to the recent progress in understanding 
the X-ray spectral characteristics of Seyfert 2 galaxies, our 
knowledge of their X-ray variability properties remains lim- 
ited. According to the standard unification scenario, the 
hard X-ray continuum should vary with large amplitude 
in a similar way to that observed in Seyfert 1 galaxies 
(Mushotzky, Done & Pounds 1993). However, in Seyfert-2s 
a large fraction of the emission must come from reprocessed 
radiation. Since a large fraction of the Fe-K line, the reflec- 
tion emission and the soft excess are likely to originate from 
regions of parsec scale-size, it follows that significantly less 
variability might be expected in Seyfert 2 objects, at least 
in those parts of the spectrum where the reprocessing makes 
a substantial contribution to the overall flux. Time variabil- 
ity studies can provide strong constraints on the geometry 
of the nucleus and the surounding region. For example, the 
time lag between the variation of the power-law continuum 
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and the Fe line flux or the reflection component would give 
information on the location and the size of the reflecting 
material. In particular if a large amount of the Fe line origi- 
nates from the accretion disk, instead of the torus, we would 
expect variations of the line flux in timescales of days. 

Here, we present the results from several RXTE archival 
observations on four well known Seyfert-2 galaxies: ESO103- 
G35, IC5063, NGC4507 and NGC7172. In another paper 
(Georgantopoulos et al. 1999) we presented monitoring ob- 
servations of Markarian 3. The extended energy range of 
the RXTE PCA detectors and the large number of observa- 
tions for each object give us the opportunity to investigate in 
detail the properties of the spectral components present in 
these galaxies such as the intrinsic power-law and the reflec- 
tion component. Our main objective is to use the variability 
exhibited in the 3-24 keV band to place constraints on the 
geometry of the nucleus and its circumnuclear matter. 



2 THE SAMPLE 
2.1 ESO 103-G35 

The HEAO-1 hard X-ray source 1H 1832-653 has been iden- 
tified with the galaxy ESO 103-G35 at a redshift z=0.013. 
Optical spectroscopy revealed a high excitation forbid- 
den line spectrum with weak broad emission line wings 
and therefore classified this object as a Seyfert-1.9 galaxy 
(Phillips et al. 1979). EXOSAT obtained the first X-ray spec- 
trum showing a power-law spectrum absorbed by a column 
of N H ~ 10 23 cm" 2 (Warwick, Pounds & Turner 1988). The 
EXOSAT observations also showed a variation in the col- 
umn density over a period of 90 days. Smith & Done (1996) 
present Ginga observations detecting an Fe line around 6.4 
keV with an equivalent width of ~ 350 eV. The spectrum 
of the source could be characterised by a heavily absorbed 
power-law continuum of slope V ~ 1.8 and Nh ~ 2 x 10 23 
cm~ 2 . Smith & Done (1996) also found evidence for the pres- 
ence of a reflection component. Turner et al. (1997) present 
the first ASCA data for this object. They resolve the Fe line 
into possibly three components (6.4, 6.68 and 6.96 keV). 
Recently, Forster, Leighly & Kay (1999) presented three 
ASCA observations separated by ~ 2 years. They find some 
marginal evidence for a double Fe line in their first observa- 
tion. Their best fit F and Nh values agree with the Ginga 
results. They also find that the power-law continuum flux 
increased by a factor of two while the Fe line equivalent 
width has decreased by about a similar amount. Two ob- 
servations of ESO103-G35 with BeppoSAX separated by a 
year, from ~2 to 60 keV (Akylas et al. 2000), show a vari- 
ation of the power-law flux by a factor of two, no variation 
of the Fe line flux while there is no strong evidence for the 
presence of a reflection component. Comparison of the his- 
torical data (Poletta et al. 1996) from HEAO-1 up to ASCA 
show a variation of the flux of a factor of about four. 



2.2 IC5063 

IC5063 is an SO galaxy (z=0.011) presenting a typical 
Seyfert-2 galaxy spectrum (Colina, Sparks & Macchetto 
1991). Scattered broad H a emission has been detected by 
Inglis et al. (1993). HSTNICMOS observations show a very 



red unresolved point source (Kulkarni et al. 1998). IC5063 
has a radio power at 1.4 GHz two orders of magnitude 
greater than that of typical Seyfert galaxies and thus it can 
be classified as a Narrow-Line Radio galaxy as well (Ulvestad 
& Wilson 1984). A Ginga observation showed a power-law 
photon index of V = 1.5 and N H ~ 2.5 x 10 23 cm -2 (Smith 
& Done 1996). Turner et al. (1997) presented two ASCA 
observations of IC5063. The spectrum is in good agreement 
with the Ginga observations. IC5063 shows short-term X-ray 
flux variability at the 90 per cent confidence level (Turner 
et al. 1997). However, the Ginga and the ASCA observation 
obtained about four years later give similar fluxes. 

2.3 NGC4507 

NGC4507 is a nearby (z=0.012) barred spiral galaxy. Optical 
spectra present high excitation, narrow emission lines clas- 
sifying it as a Seyfert-2 galaxy (Durret & Bergeron 1986). 
Ginga observations showed a flat X-ray spectrum (r ~ 1.3), 
with Nh ~ 4 x 10 23 cm~ 2 and a strong Fe line (equivalent 
width, EW ~ 800 eV) (Smith & Done 1996). Instead, OSSE 
observations (Bassani et al. 1995) showed a steeper photon 
index (r ~ 2.1 ±0.3) in the 50-200 keV energy range suggest- 
ing the presence of a strong reflection component. ASCA ob- 
servations (Turner et al. 1997, Comastri et al. 1998) showed 
a flat power-law index, with a possible variation of the iron 
line intensity and the absorption column density. The lim- 
ited spectral bandpass of ASCA did not allow to constrain 
the properties of a reflection component. Comparison be- 
tween the EXOSAT and ASCA observations show a flux 
variability of a factor of about two. 

2.4 NGC7172 

NGC7172 is an SO galaxy (z=0.0087) belonging to a compact 
group of galaxies (Hickson90). It is classified as a Seyfert-2 
galaxy on the basis of its optical spectrum (Sharpies et al. 
1984). EXOSAT obtained the first X-ray spectrum of the 
source yielding a photon index of V ~ 1.8 (Turner &: Pounds 
1989). The same photon index value was also found in a 
Ginga observation of the source (Smith & Done 1996). The 
Ginga data analysis yielded a column density of ~ 1 x 10 23 
cm -2 , and a rather weak Fe line (EW ~ 50 eV). Turner 
et al. (1997) find an EW of ~ 80 eV. The ASCA data as 
well as combination with GRO/OSSE observations (Ryde 
et al. 1997) show a flat spectrum of T ~ 1.5 and therefore 
provided evidence for spectral variability between the Ginga 
and ASCA epoch. The 2-10 keV flux has been fairly constant 
during the 1977-1989 period at a level of about 3-4 x 10~ n 
erg cm~ 2 s _1 . However, Guainazzi et al. (1998) presented 
evidence for significant short term (hours) and long term 
(months) variability using two ASCA observations separated 
by a year. They find a flux decrease by a factor of 3 over a 
year. The Fe line flux decreased by a similar amount pro- 
viding important constraints for the size of the reflecting 
material. 

3 DATA ANALYSIS 

In this work, we present the results from the analysis of 
the PCA (Proportional Counter Array) data only. The PCA 
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consists of five collimated (1° FWHM) Xenon proportional 
counter units (PCU). The PCU are sensitive to energies be- 
tween 2 and 60 keV. However, the effective area drops very 
rapidly below 3 and above 20 keV. The energy resolution 
is 18 per cent at 6.6 keV (Glasser, Odell & Seufert 1994). 
The collecting area of each PCU is 1300 cm 2 . We use only 
3 PCUs (0 to 2); the data from the other two PCU were 
discarded as these detectors were turned off on some occa- 
sions. We extracted PCU light curves and spectra from only 
the top Xenon layer in order to maximize the signal-to-noise 
ratio. The data were selected using standard screening cri- 
teria: we exclude data taken at an Earth elevation angle of 
less than 10°, pointing offset larger than 0.02°, during and 
30 minutes after the satellite passage through the South At- 
lantic Anomaly (SAA), and when the particle counts elec- 
tron0,l,2 are higher than 0.1. 

We use the PCABACKEST v2 routine of FTOOLS v 4.1.1 
to generate the background models which take into account 
both the cosmic and the internal background. The internal 
background is estimated by matching the conditions of the 
observations with those in various model files. Most of the 
internal background is correlated with the L7 rate, the sum 
of 7 rates from pairs of adjacent anode wires. However, there 
is a residual background component correlated with recent 
passages from the SAA. Therefore, the use of a second, acti- 
vation component is also employed. Comparison with blank 
fields observations shows that the above background model 
represents the real background accurately (within 1.5 per 
cent in the 3-20 keV band). 

Due to the large field-of-view of the PCA detector, 
contamination of the observed spectrum by nearby sources 
is likely. We have therefore checked the ASCA CIS 2-10 
keV images of our sources. These show no nearby sources 
(the CIS field-of-view is 40 arcmin diameter). However, 
Turner et al. (1997) point out that even the ASCA spec- 
trum of NGC4507 may be somewhat contaminated by an 
adjacent point source (present in the soft ROSAT PSPC 
0.1-2 keV band). This source cannot be resolved within the 
limited ASCA CIS spatial resolution (about 3 arcmin half- 
power-diameter) . Possible contamination at radii larger than 
the ASCA CIS field-of-view can be investigated using the 
ROSAT PSPC images of our sources. Indeed, in NGC7172 
and ESO103-G35 we find bright sources at distances from 30 
to 60 arcmin from our targets. Inspection of the NASA ex- 
tragalactic database (NED) reveals no coincidences. There- 
fore, these are most probably Galactic stellar sources and 
thus with soft X-ray spectra whose contamination in the 
RXTE band is expected to be small. 



4 THE OBSERVATIONS 

In total 43 observations with exposure time of 151.0, 46.3, 
118.6, and 78.2 ksec have been obtained with RXTE for 
ESO103-G35, IC5063, NGC4507 and NGC7172 respectively. 
The period of observations ranges from about 7 days for 
NGC7172, and 15 days for NGC4507, up to about 5 months 
for IC5063 and 7 months for ESO103-G35. The observation 
dates for each dataset together with the exposure time as 
well as the observed background-subtracted count rate in 
the full 2-60 keV PCA energy band are given in Tables 1 to 



RXTE observations of Seyfert-2 galaxies 

Table 1. The ESO103-G35 observations 



Date 


Exp (kscc) 


Ctr 


ll-Apr-1997 


14.4 


15.38±0.11 


12-Apr-1997 


10.7 


14.94±0.11 


12-Apr-1997 


15.7 


14.68±0.10 


13-Apr-1997 


6.3 


15.76±0.15 


18-Jul-1997 


6.2 


12.08±0.16 


20-Jul-1997 


8.6 


13.10±0.14 


01 T.il 1 (107 


11. 1 


1/11 O-Ln 1 1 


22-Jul-1997 


11.4 


15.92±0.12 


23-Jul-1997 


10.4 


15.36±0.11 


24-Jul-1997 


11.5 


15.32±0.11 


26-Jul-1997 


12.4 


14.30±0.11 


13-Nov-1997 


1.9 


13.69±0.27 


14-Nov-1997 


8.7 


12.84±0.13 


14-Nov-1997 


19.4 


13.20±0.09 


14-Nov-1997 


1.4 


16.10±0.32 



Table 2. The IC5063 observations 



Date 


Exp (kscc) 


Ctr 


22-Jul-1996 


18.0 


6.34±0.07 


Ol-Sep-1996 


9.8 


4.45±0.10 


Ol-Sep-1996 


1.4 


3.42± 0.24 


19-Dec-1996 


11.1 


4.61±0.09 


20-Dec-1996 


6.0 


4.73±0.12 



Table 3. The NGC4507 observations 



Date 


Exp (ksec) 


Ctr 


24-Feb-1996 


12.6 


7.26±0.08 


24-Feb-1996 


3.3 


6.96±0.16 


26-Feb-1996 


12.7 


8.33±0.08 


26-Feb-1996 


14.0 


6.99±0.07 


27-Feb-1996 


3.6 


7.46±0.15 


28-Feb-1996 


14.5 


6.75±0.08 


28-Feb-1996 


15.2 


5.39±0.09 


29-Feb-1996 


0.6 


5.83±0.35 


29-Feb-1996 


3.2 


6.84±0.07 


29-Feb-1996 


3.2 


6.73±0.16 


03-Mar-1996 


12.5 


7.70±0.08 


07-Mar-1996 


16.3 


9.58±0.07 


07-Mar-1996 


0.8 


9.04±0.32 


10-Mar-1996 


6.1 


7.69±0.11 



Table 4. The NGC7172 observations 



Date 


Exp (kscc) 


Ctr 


18-Dec-1996 


12.8 


7.32±0.09 


18-Dec-1996 


4.4 


5.42±0.14 


18-Dec-1996 


5.5 


7.61±0.12 


19-Dcc-1996 


10.0 


8.13±0.10 


19-Dec-1996 


3.2 


8.38±0.17 


23-Dcc-1996 


14.3 


6.12±0.08 


23-Dec-1996 


9.5 


4.84±0.09 


25-Dcc-1996 


11.7 


2.25± 0.09 


25-Dec-1996 


6.6 


2.47± 0.11 
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In Fig. |l| we present the 2-10 keV light curves for all 
sources using a bin size of 25ksec (time in this figure is mea- 
sured from the beginning of the first observation of each 
source). The light curves show variations on all sampled 
timescales (days/months). For example, ESO103-G35 shows 
a ~ 30 per cent variation within a few days, while IC5063 
shows a decline by a factor of ~ 2 within 4 months. The 
maximum amplitude of the observed variability is ~ 3 in 
the case of NGC7172, ~ 2 for NGC4507 and IC5063, and 
~ 60 per cent in the case of ESO103-G35 light curve. Conse- 
quently, on timescales larger than ~ 1 day, these four Seyfert 
2 galaxies show 2 — 10 keV flux variability similar to that 
observed in Seyfert 1. In fact, the light curves show signifi- 
cant variations on time scales as short as ~ a few hundred 
seconds as well (these results will be presented elsewhere). 
This is not an unexpected result, since for these Seyfert 2 
galaxies we believe that we can detect the nuclear emission 
above a few keV directly. 

4.1 Hardness Ratios 

Next, we investigate if there is any spectral variability dur- 
ing the flux variations. For this reason we extracted light 
curves at the following energy bands: 2 — 5, 5 — 7, 7—10 
and 10 — 20 keV, using a bin size of 5400 sec (this is roughly 
equal to the orbital period of RXTE ). Variations of the 
emission/absorption components in the spectrum of Seyfert 
2 galaxies are expected to contribute in a different way in 
these energy bands. For example, any absorbing column den- 
sity changes around 1 x 10 23 cm~ 2 should affect mainly the 
2 — 5 keV band, while the Fe-K emission line and the re- 
flection component variations should affect the 5 — 7 keV 
and 10 — 20 keV bands respectively. On the other hand, 
the 7—10 keV band light curve should be representative 
of the power-law continuum variations only. Using these 
light curves we calculated three hardness ratios: HRl = 

CR7-10keV /CR,2-5keV , H R2 = C Rl0~20keV / C Rr-10keV , 

and HR3 = C Rj-iokev / C R^^kev , where CRe\-e2 is the 
count rate of the El — E2 band light curve. Under the hy- 
pothesis that there are no spectral changes, the HRs should 
remain constant. 

Fig. ^| shows the HR plot for the sources as a function of 
the normalised C 'R7-iokeV (the normalisation is necessary 
in order to compare the sources since they have different 
luminosities). NGC7172 shows the largest amplitude vari- 
ations in the 7—10 keV band. Within the common range 
of the normalised CRr-iokeV values, the HRs have similar- 
values for all objects except for NGC4507 (open diamonds 
in Fig. ^). Both HR2 and HR3 for this source have a large 
value. In fact we could not calculate HR1 as its 2 — 5 keV 
light curve has a mean value of almost zero. The similarity 
of the HR values for ESO103-G35, IC5063 and NGC7172 
suggests that they have similar spectra, while the NGC4507 
spectrum appears to be much "harder". This is probably 
due to the fact that NGC4507 has a larger absorbing col- 
umn density, which reduces mainly the count rate in the 
"soft" bands, and hence results in larger HR values. 

In order to investigate the presence of spectral varia- 
tions we fitted the data for each galaxy with a model of the 
form: log(HR) = a + b* log(CRr-iokev) ( m fitting the data 
we took account of the errors in both HR and C7?7_iofc e y )• 
In all cases, a line in the log-log space fits the data well (ie 



the best fit \ 2 values can be accepted with confidence larger 
than 5 per cent). If the HR values remain constant during 
the flux variations, we expect b (ie the slope of the line) to 
be consistent with zero. We discuss below the best fit results 
for each galaxy. 

ESO103-G5 shows no HRl variations (b H m = 0.05 ± 
0.06). On the other hand, both HR2 and HR3 are variable: 
b HR 2 = -0.23 ± 0.06 and b HR3 = -0.10 ± 0.04. The HR2 
variations suggest that the spectrum of the source above 7 
keV steepens as the source brightens. The HR3 ratio shows 
the opposite behaviour. It implies a "flattening" of the 5 — 
10 keV spectrum as the source brightens. In the case of 
IC5063, both HRl and HR2 are variable (b H m = -1.0 ± 
0.2, and bHR2 = —0.40 ± 0.10) while HR3 remains constant 
(b Hm = 0.06 ±0.07). Both HRl and HR2 indicate that the 
source's spectrum becomes softer as it brightens. NGC4507 
shows a different behaviour. HR2 remains constant (bnR2 

0,04 ± 0.03), while HR3 is significantly variable (bnR3 = 
0.45 ±0.08). The spectrum below 10 keV appears to become 
softer as the source flux increases. NGC7172 shows a similar 
behaviour. HR2 remains constant (pHR2 = —0.04 ± 0.03), 
HR3 shows marginal changes (bHR3 = 0.06 ± 0.03) while 
HRl shows more significant variations (pHR\ = — 0.12±0.4). 
These variations indicate again a steepening of the spectrum 
below 10 keV as the source flux increases. 



5 SPECTRAL ANALYSIS 

In our spectral analysis, we used only data between 3 and 20 
keV, where the effective area of PCA, and hence the signal- 
to- noise ratio, is highest. By discarding the data below 3 keV 
we also avoid the complications associated with the soft X- 
ray excess in these sources (see Turner et al. 1997). The PCA 
data from each observation were binned to give a minimum 
of 20 counts per channel (source plus background). The spec- 
tral fitting analysis was carried out using the XSPEC v. 10 
software package on the basis of "joint fits" to our RXTE 
observations. All errors correspond to the 90 per cent confi- 
dence level for one interesting parameter. All energies quoted 
refer to the emitter's rest-frame. 

We have employed two group of models (see below) in 
order to characterize the spectra of the objects at each obser- 
vation. Within each group, in order to assess the significance 
of new parameters added to the initial model of the group 
we have adopted the 95 per cent confidence level using the 
F-test for additional terms. To compare models with dif- 
ferent number of parameters within each group and models 
between the two groups we have used the "likelihood ratio" 
(Mushotzky 1982) and accepted values of the ratio larger 
than 100 as showing significant improvement in the model 
fits. 

Note that the RXTE background estimation software 
calculates errors on the simulated background spectrum as- 
suming Poisson statistics. This is an overestimate of the 
statistical error on the background spectrum, which is es- 
timated using a large amount of data. In fact, for faint ob- 
jects (such as AGN) the error that the standard software 
calculates for the background is comparable to the error in 
the total observed rate. In practice, the background spec- 
trum errors should be negligible compared to the errors on 
the observed spectrum. Since the XSPEC software estimates 
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Figure 1. Light curves from the PCA, 0,1,2 PCUs, 2-10 keV band for the four sources. Bin size is 25ksec, and the time axis in seconds 
from the beginning of the first observation for each source. 
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Figure 2. Hardness ratios vs normalised 7-10 keV count rate for the four sources (for the definition of HR1, HR2 and HR3 see section 
4.1). The hardness ratios were estimated using 5400 sec binned light curves. Open circles are for ESO103-G35, filled triangles for IC5063, 
open diamonds for NGC4507 and filled squares for NGC7172. 



the error in the net spectrum (source - background) by com- 
binning the errors in the observed spectrum with those of 
the simulated background, the reduced \ 2 values in many 
cases are very small. However, the F and "likelihood ratio" 
tests are not affected by this. Finally, due to the large er- 
rors, the confidence regions that we estimate for the best fit 
parameter values are almost certainly enlarged. 



5.1 Power Law and Gaussian Line Models. 

Following previous Ginga and ASCA results, we first employ 
a simple spectral model consisting of a power-law contin- 
uum, with photon spectral index F, modified by absorption 
in a column density, Nh, of cool neutral material. A Gaus- 
sian line was also included to account for Fe-K emission. 
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The line energy and the width were allowed to vary freely. 
However, when the width was found to be smaller than the 
RXTE spectral resolution, we fixed it at au ne = 0.1 keV. 
The values of the Nh, the photon index and the normaliza- 
tion of the Fe line are free parameters but tied to the same 
value for all observations. The normalization of the power 
law is allowed to vary freely. The results of fitting this spec- 
tral model are presented as model A in Tables 5, 6, 7 and 
8. Next, we allow for the column density and the power-law 
index to vary freely (models B and C in the same tables). 
Finally, in model D we allow for the normalization of the 
Gaussian line to vary freely while now both the power-law 
index and the column density are tied to the same value be- 
tween all the observations. We discuss the results for each 
galaxy in turn. 

ESO103-G35. It shows a large column density (Nh ~ 
1.4 x 10 23 cm -2 ) in excellent agreement with the Ginga 
observations. However, the spectral index is flat (r ~ 1.6, 
model A) in contrast with the Ginga results. Moreover, the 
best fit energy line is very low (E ~ 6 keV) and marginally 
resolved, ie broader than the energy resolution of PCA, hav- 
ing a a of ~ 0.7 keV. The EW of the line does not remain 
constant. It ranges from 490lgg to 661 eV. The data show 
strong evidence for spectral variability, since models B,C and 
D give a significantly better fit to the data when compared 
to model A. When the line normalisation is untied (model 
D), we obtain a A% 2 = 68 for 14 addtitional parameters. 
A larger A\ 2 value is obtained when the column density is 
untied and an even larger A% 2 when the power law index 
is untied (A^ 2 = 107 and 162 respectively for 14 additional 
parameters). Based on the likelihood ratios, model C gives 
the best fit to the data. 

Following Forster et al. (1997), we tried an additional 
model. We added an absorption edge component (model 
EDGE in XSPEC) to model A and kept the Fe line energy 
fixed at 6.4 keV. We obtain a highly statistically significant 
improvement in x 2 compared to model A (A^ 2 ~ 143 for 
one additional parameter). The photon index becomes now 
steeper (r = 1.68±g;g| and N H = 1.71° j x 10 23 cm -2 ). The 
best-fit energy for the edge is 6.8 jig io keV and the optical 
depth is r = 0.2 ± 0.02. 

IC5063. Model A does not provide a good fit in this 
case. The column density is 2 x 10 23 cm _2 while the power- 
law is flat (r ~ 1.6) both in good agreement with the Ginga 
results. The EW of the Fe line varies between 298±H and 
462±^j eV. As before with ESO 103-G35, we find highly 
significant evidence for spectral variability. Models B and C 
give a significant % 2 improvement compared to model A (in 
this case model D does not improve the fit significantly) . The 
best fit is obtained when the photon index is untied (model 
C): r = 1.15— 1.70. However, the Fe line energy is low (E ~ 
6.12 keV with a ~ 0.5). We therefore included an absorption 
edge in model C. The line energy is marginally improved, 
E = 6.19±g;J| keV, with a = 0.48. The edge energy was 
fixed at 6.9 keV while the optical depth found was r = 0.07. 
However, the inclusion of the edge did not improve the fit at 
a statistically significant level (x 2 = 175.2/235, as compared 
to x 2 = 177.2/236 for model C). 

NGC4507. Model A results in a high absorbing column 
(~ 4.5 x 10 23 cm -2 ) and a flat spectrum (T ~ 1.4), in agree- 
ment with the Ginga observations. The line energy is ~ 6.16 
keV with the line width being unresolved (a ~ 0.33). The 



EW ranges from 353±|i up to 655±|f eV. All models B,C 
and D give a significant \ 2 improvement when compared to 
model A, indicating again the presence of significant spec- 
tral variability. Based on the likelihood ratios, it is model B 
this time that provides the best fit to the data. 

NGC7172. Model A shows a large column density 
(Nh ~ 1-1 x 10 23 cm -2 ) in agreement with previous Ginga 
results. The power-law has a value of F ~ 1.8, also in good 
agreement with the Ginga results, and typical of the intrin- 
sic AGN spectral index. The Fe line has an energy of ~ 6.24 
keV. The line width is consistent with the line being narrow. 
The resulting EW varies from 96±22 to 360± 28 eV in the nine 
observations. Model D does not improve the fit, while mod- 
els B and C give a statistically significant improvement to 
the x 2 value of the fit. When the column density becomes a 
free parameter (model B) we obtain a A^ 2 ~ 30 for 12 addi- 
tional parameters. When the photon index becomes untied 
(model C), we obtain an equally good x 2 with T varying 
between 1.7 and 2.0. This suggests that the two parameters 
(photon index and column density) are degenerate and it is 
therefore difficult to disentangle the origin of the observed 
spectral variation. We plot the variations of column density 
and of the photon index as a function of the 3 — 20 keV flux 
in Fig. ^ and ^ respectively. We can see that the changes in 
the spectrum of this source can be explained by a decrease of 
Nh by ~ 30 per cent as the source flux increases by a factor 
of ~ 3.5. Alternatively, a F increase by ~ 0.3 (from ~ 1.65 
up to ~ 1.95) can explain the spectral changes equally well. 

5.2 Power Law, Gaussian Line and Reflection 
Component Models. 

Although the power-law plus gaussian line prescription de- 
fined above gives an acceptable fit in terms of the x 2 statis- 
tic, there is evidence from the flat power-laws derived for 
NGC4507, IC5063 and ESO103-G35 for the presence of a 
strong Compton reflection component. The next step in the 
current analysis was therefore to include a reflection com- 
ponent in the spectral modeling. Specifically we use the 
PEXRAV model (Magdziarz & Zdziarski 1995) in XSPEC. 
This model calculates the expected X-ray spectrum when a 
source of X-rays is incident on optically thick, neutral (ex- 
cept hydrogen and helium) material. We assume that both 
the reflection component and the power-law are absorbed 
by the cold gas column density. The strength of the reflec- 
tion component is governed by the parameter R, represent- 
ing the strength of the reflected signal relative to the level 
of the incident power-law continuum. We set R=l (which 
corresponds to a 2n solid angle subtended by the optically 
thick material). We fix the inclination angle for the disk at 
i = 60°, since the shape of the reflection spectrum below 
20 keV is relatively independent of the inclination angle. 
We also fix the energy of the exponential cutoff at 300 keV. 
Both the iron and light element abundances were kept fixed 
at the solar abundance values. 

We added this reflection component to all models A,B,C 
and D, creating a new set of models (A',B',C and D'). In all 
new models, the normalisation of the reflection component 
is tied to a single value across the set of observations. We 
also set the spectral index of the power-law continuum to 
r = 1.9 (in all models apart from model C). As before, in 
model A' the power-law normalisation is allowed to vary in 
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each observation while the column density and the Fe line 
normalization are tied between the different observations. In 
model B' we untie the column density, in model C we untie 
the photon index and in model D' we let the gaussian line 
normalisation free. Finally, we let the reflection component 
normalisation free (model E'). The results are presented in 
tables 9 to 12. 

ESO103-G35. In the case of ESO103-G35, the inclu- 
sion of a reflection component gives a statistically signifi- 
cant improvement to the data fit. For example, we obtain a 
Ax 2 = 163 (with the addition of one extra parameter) when 
we include the reflection component in model A. This is in 
agreement with the findings of Smith & Done (1996) using 
Ginga data. Furthermore, all models (B',C',D' and E') sig- 
nificantly improve the fit when compared to model A'. The 
best fit to the data is obtained when the photon-index is al- 
lowed to vary with time (model C, T ~ 1.8 — 2). An equally 
good fit is also obtained when the reflection component nor- 
malisation is allowed to vary with time (model E'). Both 
these models give a better fit than model C, showing again 
that the inclusion of the reflection component significantly 
improves the fit. 

In Fig. |§] we present the variations of the power-law and 
the reflection flux between the different observations (model 
E'; time in this figure and in Fig. ^ as well is measured 
in seconds from the mission launch date). We see that the 
power-law component shows rapid flux variations (up to a 
factor of 50 per cent) within days. Similar variations, but not 
well correlated with the power-law variations, are observed 
in the reflection component flux as well. In Fig. ^ we plot 
model's C best fit V values as a function of the 3 — 20 keV 
flux. A change of ~ 0.15 (from V ~ 1.8 up to F ~ 2) can 
explain the spectral variability in this case. 

However, the best fit line energy in models C and E' is 
abnormally low (E ~ 6.15 keV). We therefore introduced an 
absorption edge to model E'. The resulting \ 2 ls 505/714, 
which indicates that, compared to model E', the improve- 
ment to the fit is statistically significant. The best fit line 
energy now becomes 6.23±0.10 keV, roughly consistent with 
the line energy of neutral Fe. The best fit edge energy is 
E e dge = 7.69 ±0.14 keV (ionised Fe) while the optical depth 
is t = 0.1 ± 0.01. Similar results are obtained when we add 
a warm edge component to model C 

IC5063. In IC5063 the inclusion of a reflection compo- 
nent does not represent a statistically significant improve- 
ment. Models C, E' and C (which does not include a reflec- 
tion component) give an equally good fit to the data. How- 
ever, the inclusion of the reflection component (for example 
model E') yields a line with an energy of E = 6.25 ± 0.10 
keV, which is roughly consistent with the energy of line 
from neutral Fe. In Fig. ^ the variations of the power-law 
and the reflection flux are presented (model E'). The power- 
law component flux decreases by a factor of two with time. 
In contrast, the reflection signal, although variable as well, 
does not correlate with the power-law flux. Alternatively, a 
change of F between roughly 1.6 and 2.2 (model C, see Fig. 
^|) yields a similarly good fit to the data, and results in the 
same best fit line energy value. However, the photon index 
does not correlate well with the continuum flux. 

NGC4507. In this case, the use of a reflection compo- 
nent improves the fit dramatically. Models B' and C give a 
much better fit to the data than model B. These two models 
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Figure 3. ESO103-G35: The 3-20 keV flux of the power-law (up- 
per panel) and of the reflection component (lower panel) as a 
function of time from the mission's launch (90 per cent errors are 
plotted) 



(untied column densities or photon indices) fit equally well 
the data: in the first case the Nh varies roughly between 
4.8 and 5.3 x 10 23 cm -2 , while in the second case F varies 
between 1.9 and 2.0. 

The line energy is again low (E ~ 6.2 keV) and therefore 
we included an absorption edge component in both models 
B' and C above. We obtain a statistically significant im- 
provement in the fit in both cases. However, the best fit is 
now clearly obtained in the case of a column density varia- 
tion (x 2 =643/794 as compared to x 2 =655/794 for an untied 
r). The edge energy was frozen at 6.9 keV while its optical 
depth was found to be 0.16 ±0.02. The line energy was kept 
fixed at 6.4 keV (<r=0.1 keV). In Fig. |H] we plot the varia- 
tions of Nh as a function of the 3-20 keV flux (model B' with 
edge). A ~ 10 per cent decrease in the absorpting column 
while the source flux increases by a factor of 2 can explain 
the observed spectral changes. 

NGC7172. In this case, the inclusion of a reflection 
component does not improve the fit at a statistically sig- 
nificant level. Models B' and C present the best fit to the 
data, however the goodness of fit is comparable to that of 
the models B and C. The lack of a reflection component is 
rather expected given the steep power-law index found ear- 
lier in the case of the simple power-law and a Gaussian line 
model. Note that the Ginga results required the presence of 
a reflection component only marginally. 



5.3 Soft excess models 

The X-ray spectra of Seyfert-2 galaxies usually show some 
level of soft X-ray emission despite the strong photoelectric 
absorption (eg Turner et al. 1997). This soft excess could 
be due to either a strong star-forming component or scat- 
tered emission from the nucleus. The contribution of the 
soft excess in our case (ie above 3 keV) should be negligi- 
ble. Nevertheless, we tried to investigate whether a constant 
soft excess combined with a power-law component with vary- 
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Figure 5. The variation of the column density in NGC4507 (model B' with edge) and NGC7172 (model B) as a function of the 3-20 
keV flux 
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Figure 6. The variation of the photon index in ESO103-G35 (model C with edge), IC5063 (model C) and NGC7172 (model C) as a 
function of the 3-20 keV flux 



ing normalization may explain the strong spectral variabil- 
ity observed in our data. Therefore, we include a power-law 
component with the same photon index as the hard X-ray 
power-law and with a free normalization (but tied between 
all the observations) in the models above. 

In NGC7172 we add the soft power-law to model A 
since a reflection component is not needed by the data. We 
obtain x 2 =303.8 ie identical to that of model A, despite the 
addition of one extra parameter. We conclude that an extra 
soft component is not needed by the data. In the other three 
cases we add the soft power-law component to model A'. 
In the case of ESO103-G35 we obtain X 2 =568.5/729 which 



represents a statistically significant improvement to model 
A'. However, the fit is worse when compared to models C 
and E'. In NGC4507 we obtain \ 2 = 738/806 similar to 
model A' with the soft power-law normalization being zero. 
Finally, in the case of IC5063 we obtain an identical \ 2 to 
model A' (269/239) while the soft power-law normalization 
is again zero. The above confirm that the presence of any 
soft excesses do not affect our data above 3 keV, and they 
cannot account for the observed spectral variability. 
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Figure 4. IC5063: the 3-20 keV flux of the power-law (upper 
panel) and of the reflection component (lower panel) as a function 
of time from the mission's launch (90 per cent errors are plotted). 



6 SUMMARY 

ESO103-G35 presents a complex spectrum. Models C'and 
E', with the addition of a warm edge, give the best fit to 
our data. Our best fit edge energy is 7.69 ± 0.14 keV with 
r = 0.1 ± 0.01. The best fit edge energy matches the re- 
sults from the 1988 and 1991 Ginga observation (Warwick 
et al 1993, Smith and Done 1996), and also the 1996 ASCA 
observation (Forster et al 1999) . The energy of the edge im- 
plies the presence a significant amount of ionised iron. The 
optical depth of the iron edge was observed to increase be- 
tween the 1994 and 1996 ASCA observations (Forster et al 
1999), and the edge line in the former observation was con- 
sistent with that from neutral iron. Our best fit edge optical 
depth is significantly smaller than the respective ASCA re- 
sult (Forster et al. 1999). However, the quality of the present 
data set is not good enough to investigate if there exist any 
changes in the optical depth and/or the energy of the edge 
during the RXTE observations. The spectral variability can 
be explained by a variable spectral index (model C, Fig. |^). 
Although these variations are subtle, they are highly sta- 
tistically significant. The photon index is correlated with 
the 3-20 keV flux in the sense that the spectrum steepens 
with increasing flux. The same variations can be equally 
well explained by a variable reflection component (model 
E', Fig. ||). Interestingly, Akylas et al. (2000) find no reflec- 
tion component in the BeppoSAX observations despite the 
large passband (~2-60 keV) and the good signal-to-noise ra- 
tio of their observations. This may suggest that the relative 
normalization of the power-law and the reflection compo- 
nent must change in different epochs, and supports the in- 
terpretation that the spectral variations seen by RXTE are 
caused by a variable reflection component. Furthermore, in 
this case we can also explain naturally the fact that although 
the HR2 ratio of ESO103-G35 is variable, the HRl ratio re- 
mains constant. 

IC5063 can be well fit with a simple spectrum consist- 
ing of a power-law and a Gaussian line but the energy of 
the line is abnormally low (~ 6.12 keV). The addition of 



an edge does not yield a reasonable value for the energy of 
the Fe line. However, the addition of a reflection component 
(models C and E'), although it does not improve the X* , 
yields an Fe line with an energy of 6.25 ± 0.1 keV compati- 
ble within the errors with the energy of cold Fe line. IC5063 
shows strong evidence for spectral variability (Fig. ^). The 
spectral variability can be explained either by the reflection 
component being variable between the observations (model 
E') or by the photon index being variable (model C, T ~ 1.6 
- 2.2). According to model E' best fit results there is no cor- 
relation between the reflection component and the intrinsic 
power law flux variability (Fig. [|). Similarly, in the model C 
case, there is no clear correlation between the photon index 
and the continuum flux (Fig. fy. 

NGC4507 presents a similar spectrum to ESO103- 
G35. We find strong evidence for a reflection component 
(note that such a component was not evident in neither the 
Ginga nor the ASCA data). We also find evidence for a cold 
Fe edge (with an energy fixed at 6.9 keV), in excess of the ab- 
sorption edge associated with the reflection component. As- 
suming the iron cross section given by Leahy and Creighton 
(1993), the optical depth of the edge (r = 0.16±0.02) implies 
an equivalent hydrogen column density of ~ 2 x 10 23 cm" 2 . 
This is smaller than our best fit Nh value of 5 ± 1.0 x 10 23 . 
The edge component is not required by the data, but with 
its addition we can get a good fit to the data with the iron 
line energy fixed at 6.4 keV. However, there is a possibility 
that residuals in the PCA calibration matrix around 5 keV 
due to Xe may still be present. In such a case it is possi- 
ble that the line energy obtains somewhat low values due to 
these residuals. Therefore, we believe that the detection of a 
cold Fe edge in this object remains uncertain. We find strong 
evidence for spectral variability (Fig. |2|). The best model fit 
to the data (model B') suggests flux correlated column den- 
sity changes: the column becomes lower with increasing flux 
(Fig. ||). This model can explain naturally the HR3 varia- 
tions observed for this object, and at the same time the lack 
of HR2 variations. 

NGC7172 is the only object which shows no evidence 
for a reflection component. Ginga found marginal evidence 
for such a component while in BeppoSAX data a reflection 
component is needed only if the spectral index is fixed to the 
canonical F — 1.9 value (Akylas et al. 2000). Similarly, if F 
is fixed to 1.9, models with reflection (B' and C') can fit the 
present data set well but simple power law models (B and 
C) can fit the data equally well. There is no evidence for an 
edge, and spectral variability is observed again. The spectral 
fitting cannot discriminate between variations in the photon 
index or the column density: as the flux increases, either the 
column drops (14 — 11 x 10 23 cm -2 ) or the photon index 
steepens (1.66 — 1.96) in a similar fashion to ESO103-G35 
and NGC4507. 



7 DISCUSSION AND CONCLUSIONS 

The RXTE observations of ESO103-G35, IC5063, NGC4507 
and NGC7172 show that, in all objects, the X-ray flux is sig- 
nificantly variable on time scales larger than ~ a day, with 
an amplitude similar to that observed in Seyfert 1 galax- 
ies. This result on its own suggests that, in these objects, 
we are seeing directly the nuclear emission. This emission 
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is transmitted through a heavily absorbing material with a 
density of ~ 1 — 5 x 10 23 cm -2 . Apart from the flux vari- 
ability, we searched for the presence of spectral variations, 
taking advantage of the fact that the present data set con- 
sists of many observations for each object all made with the 
same instrument. In the past, spectral variability studies in 
Seyfert 2 galaxies were not conclusive mainly because it is 
difficult to compare the results from different spectral fits 
performed across various instruments with different band- 
passes. We detected significant, flux related variations in 
the hardness ratios that we calculated for all objects (Fig. 
^|). This result suggests that spectral variability may be a 
ubiquitous feature in Seyfert-2 galaxies. 

In principle, the observed hardness ratio variations 
could be an artefact of poor background subtraction. How- 
ever, the HR variations do not show the same behaviour 
in all sources, and are well correlated with the source flux. 
Therefore, if these variations are due to background sub- 
traction errors, those errors should be different at different 
epochs and should operate in such a way so that the result- 
ing ratios should mimic the flux correlated variations seen in 
Fig. [| We consider this possibility unlikely to be the cause 
of the spectral variations that we observe. 

Still, the explanation for the observed spectral variabil- 
ity is model dependent and is based on the spectral fitting 
results. In some cases, these results are somehow uncertain. 
For example, in NGC7172 and NGC4507, significantly differ- 
ent Nh values are observed in just one observation, namely 
the one with the highest flux (Fig. |B|). Although the hardness 
ratios clearly show that the variability in both NGC4507 
and NGC7172 is attributed to more than one observation, 
we investigated whether the spectral fitting results could be 
affected by problems associated with the background sub- 
traction in the highest flux observations. Hence, we exam- 
ined the data below 3 keV in layer 1 where the signal is 
mostly due to particles as the effective area of the PCA is 
negligible. We also examined the signal in the full energy 
band of layers 2 and 3 which are again sensitive to particles 
only. We find that the particle spectrum was consistent with 
the background model predictions, especially below 5 keV. 
In conclusion, we believe that the spectral variability is real 
in all four objects, and that our model fit results can be used 
in order to explain its origin. 

The origin of the spectral variability appears to be dif- 
ferent in each object. We discuss below briefly the impli- 
cations of our results on the processes that operate in the 
nuclei of Seyfert 2 galaxies. 

7.1 Photon index changes 

In three cases, namely ESO103-G35, IC5063 and NGC7172, 
models with a variable photon index (model C for the first 
two sources and model C for NGC7I72) can fit the data 
well. Turner et al (1998) have found significant photon index 
variability in many Seyfert 2 galaxies by comparing histor- 
ical X-ray observations. Significant photon index variability 
has also been observed in many Type 1 objects (eg MCG 6- 
30-15, Lee et al. 1999; NGC4051, Guainazzi et al. 1996; Mrk 
766, Leighly et al. 1996; NGC4151, Warwick et al. 1996). 
In all these objects, the power law becomes steeper with 
increasing flux. The same pattern is followed by the pho- 
ton index with flux in our case as well. In ESO103-G35 and 



NGC7172 we observe F increasing by ~ 0.2 and ~ 0.3 for 
a ~ 1.6 and ~ 3 flux increase respectively. In IC5063, V in- 
creased by ~ 0.6 while the source doubled its flux, but is far 
from clear in this case that the V variations are correlated 
with the source flux. 

According to standard Comptonization models, hard X- 
rays in AGN are produced through the scattering of the 
soft UV photons from the accretion disk on a hot (> 40 
keV) corona above the accretion disk. Opacity changes in the 
corona can give rise to significant spectral variability, even 
if the total luminosity of the hard X-rays remains constant 
(Haardt et al 1997). If the corona opacity is not dominated 
by pairs, as the optical depth of the corona, r, increases the 
intrinsic spectrum becomes steeper (as long as r < 1). As 
a result, as long as F < 2, the the photon index correlates 
with the photon flux: a ~ 2.5 increase in the 2 — 10 keV 
flux can result in a Ar ~ 0.4 increase in the photon index 
(Fig. 7 in Haardt et al 1997). The best fit V values that we 
find are less than 2 in almost all cases and the AF changes 
that we detect for the observed source flux variations are in 
broad agreement with the predictions of this model. 

We note in passing that Seyfert 2 galaxies with "flat 
spectrum" (ie Smith & Done 1996) may not be intrinsically 
flat. If those sources also exhibit flux related photon index 
variations, then values as flat as ~ 1.5 — 1.6 can be expected 
according to Haardt et al. (1997). Monitoring observations 
of these sources could resolve the issue. 

7.2 Absorption column changes 

In NGC4507 and NGC7172, models with a variable absorp- 
tion column density (models B' and B respectively) give a 
good fit to the data. In fact, model B' gives the best fit to 
the NGC4507 data. In both cases, the column density is an- 
ticorrelated with the X-ray flux: a~ 5 x 10 22 cm _2 decrease 
in Nh is observed while the source flux increases by a factor 
of ~ 2 and ~ 3 respectively (Fig. |^) . Absorbing column vari- 
ations in Seyfert 2 galaxies have been observed in the past as 
well (eg NGC526A, Turner & Pounds 1989; NGC7582, War- 
wick et al. 1993; ESO103-G35, Warwick, Pounds and Turner 
1988). In the first two cases significant column changes on 
a time scale of years were detected, while two EXOSAT ob- 
servations of ESO103-G35 revealed a significant drop in the 
column density within 90 days. 

The flux correlated Nh variations can be explained if 
the obscuration is not due to a uniform torus, but due to 
individual clouds rotating around the nucleus instead. We 
may expect clouds of different column density to intercept 
the emission from the nucleus at different epochs. In this 
case, a small increase of the column density will be associ- 
ated with a decrease in the total flux. This model requires 
both high velocities for the covering clouds and a small X- 
ray source size. For example, in the case of NGC7172, the 
lowest flux observations (last two points in Fig. ji]) have sim- 
ilar Nh values and are ~ 4.5 and ~ 5.5 days apart from 
the highest flux observation (third point in Fig. |l|) . Assum- 
ing that a cloud of Nh ~ 3 — 5 x 10 22 cm _2 and velocity of 
~ 10 4 km/sec covered the X-ray source in 4.5 — 5.5 days, 
then the size of the X-ray source should be smaller than a 
few (~ 3.5 — 4.5) light hours. 

The Nh variations could also be explained by partial 
photo-ionisation of the obscuring screen. Varying degrees 
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of photo-ionisation caused by intrinsic luminosity variations 
of the source can result in column density changes. This 
model can explain succesfully the variable X-ray absorption 
in MR2251-178 (Halpern 1984, Pan et al. 1990). Recently, 
Comastri et al. (1998), presented strong evidence for the 
existence of photonionised material with a column density 
of ~ 8 x 10 22 cm" 2 in NGC4507. In order for the absorb- 
ing column to adjust itself instantaneous to the continuum 
flux variations (and hence explain the rather "continuous" 
variations shown by its HR2 plot in Fig. []), the absorbing 
material should be close to the central source and in a thin 
shell (so that its density will be larger than ~ a few 10 8 
cm -3 , Netzer 1993). Interestingly, Turner et al. (1998) have 
recently presented evidence for the existence of absorbing 
material close to the central source in Seyfert-2 galaxies. 
Neither the ASCA nor the RXTE data show evidence for 
a warm edge in the spectrum of NGC4507. For an ionised 
material with a column density of ~ a few x 10 23 cm _2 the 
respective edge should have an optical depth of less than 0.1, 
and hence difficult to detect. On the other hand, ESO103- 
G35 shows a warm iron edge in its spectrum but no in- 
dication of absorption column variations. Perhaps in this 
case, the small amplitude continuum variability together 
with possible variability of other spectral components (ie 
P, reflection component normalisation) render the detection 
of any absorption variations difficult. Finally, in the case of 
NGC7172, ASCA observations did not show any evidence of 
photo- ionised material (Guainazzi et al 1998). The 2 — 10 
keV flux of the source during the ASCA observation was 
~ 1.3xl0 -11 erg cm _2 s _1 , which is similar to the flux during 
the last RXTE observation. Their best fit F value 1.5 ±0.15, 
which is consistent with our estimate (1.66 ± 0.15, Table 8). 
This result supports the idea that in NGC7172 flux related 
r changes (ie model C) rather than Nh variations (model 
B) are responsible for the spectral variability. 

The detection of Nh variations in just one source indi- 
cates that distribution of the absorbing material in Seyfert 
2 galaxies may not be the same in all objects (eg Turner et 
al. 1998). Material away from the nucleus is not expected to 
be ionised and will not respond to the continuum variability. 
As a result, Nh in these objects will remain constant with 
time. 

7.3 Reflection component and Iron Line 

Model E' (a variable reflection component) can explain the 
spectral variability observed in ESO103-G35 and IC5063. If 
true, this would imply that the cold material responsible 
for reflection is close to the nucleus eg an accretion disk as 
suggested by Nandra et al. (1997) for Seyfert 1 galaxies. In 
this case we expect a positive correlation and no appreciable 
lags in the reflection response to the continuum flux varia- 
tions. Contrary to this, we find that the reflection variations 
are not correlated with the continuum variations, especially 
in the case of IC5063 (Fig. ^, ^). One can assume, rather 
arbitrarily, that there is after all a lag and therefore the 
reflection component is responding to past, unobserved con- 
tinuum variations. Further evidence for the existence and 
possible variability of the reflection component can be pro- 
vided by the iron line results. 

All sources show a strong, unresolved emission line at 
~ 6.4 keV. The only exception is ESO103-G35 where the line 
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Figure 7. The equivalent width of the 6.4 keV line plotted for 
all objects as a function of the 2-10 keV flux. Solid squares are 
for ESO103-G35, solid triangles for IC5063, open diamonds for 
NGC4507 and open squares for NGC7172 (90 per cent errors are 
plotted) 



may be marginally resolved(<r ~ 0.5 keV, models C and E'). 
This is consistent with the results of Turner et al. (1997) 
for this source. The intensity of the line though remains 
constant. This is indicated by the fact that model D' (a 
variable iron line) does not improve the fit when compared 
with model A' in IC5063 and NGC4507. The same holds 
for model D when compared with model A for NGC7172. 
Fig. 0, shows a plot of the line EW as a function of the 
source 2 — 10 keV flux. The equivalent widths were estimated 
using results from the best fit of model B' for NGC4507 and 
NGC7172, model C for IC5063 and the same model plus a 
warm iron edge for ESO103-G35. Since all models require 
a constant line, in all cases the EW decreases as the source 
flux increases. The errors on the EW are small enough to 
see that the EW variations are significant. 

The fact that the lines are narrow (except for ESO103- 
G35) and constant implies that they are produced by the 
obscuring torus. For IC5063, NGC4507 and NGC7172, the 
EW of the iron line can be as large as ~ 300 eV in the low- 
est flux observations (Fig. ^). The measured column densi- 
ties in the line of sight toward those sources is ~ 3 x 10 23 , 
~ 5 x 10 23 and ~lx 10 23 cm~ 2 respectively. For these N H 
values, we expect an emission line with equivalent width of 
~ 120 eV, ~ 200 eV and ~ 40 eV respectively (Ghisellini et 
al. 1994). The larger than expected EW suggests that there 
exists a source of line photons other than from the line-of- 
sight absorbing material. A 6.4 keV iron line can also be pro- 
duced by fluorescence from an acrretion disk. The expected 
angle- averaged EW in this case is ~ 110 eV. Consequently, 
in NGC4507 and IC5063 the iron line that we observed is 
probably the sum of two lines - most of the line is produced 
by the torus and there is additional contribution by the ac- 
cretion disk itself. Any variations of the accretion disk line 
would be diluted by the presence of the torus line and hence 
difficult to detect. Furthermore, the lack of line variability 
indicates that the torus should be in a distance larger than 
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approximately two light weeks and approximately six light 
months in NGC4507 and IC5063 respectively. 

For NGC7172, the discrepancy between the EW ob- 
served in the two lowest flux observations of this source 
(~ 300 eV), and the EW expected in the case of a torus 
plus an accretion disk line (~ 150 — 200 eV) remains. Fur- 
thermore, the Ginga and ASCA results (Smith & Done 1996, 
Turner et al. 1997) showed a cold iron line with EW ~ 50 
and ~ 70 eV respectively (these estimates had large errors 
associated with them). The continuum flux was ~4x 10 -11 
erg s _1 cm -2 in both cases. Despite the large errors, for this 
source flux level, the equivalent widths are in agreement with 
these expected according to EW vs flux plot for this source 
in Fig 7. Guainazzi et al. (1998) report the same EW when 
the source was ~ 4 times fainter. Their estimate appears 
to be inconsistent with our results, and suggests that the 
line follows the continuum on time scales of ~ a year. The 
RXTE observations were made within a week. One explana- 
tion could then be that the line is mainly produced by the 
torus which is located at a distance larger than ~ a light 
week and smaller than ~ a light year. 

ESO103-G35 was observed three times by ASCA 
(Forster et al. 1999). In all cases, the source flux was lower 
than during the RXTE observations. The EW of the iron 
line was higher (~ 400 — 200 eV), in agreement with the 
trend shown in Fig. 0. As in the case of NGC7172, the large 
EW observed by ASCA requires an extra source of line pho- 
tons apart from the torus and the accretion disk. Forster 
et al. (1999) suggest various possibilities in order to explain 
this. The fact that the iron line appears broad in the RXTE 
observations is a direct indication that a part of the line 
is produced by fluorescence from the accretion disk. In fact, 
ESO103-G35 is the only source were a model with a variable 
iron line (model D') improves significantly the goodness fit 
when compared to model A'. This indication of line vari- 
ability supports the idea that the spectral variability in this 
object is due to a variable reflection component (model E'). 
Such a model could also explain naturally the fact that HR2 
in the case of ESO103-G35 is variable while HR1 remains 
constant. The possible detection of line and reflection com- 
ponent variability in this object only indicates that the ori- 
entation of the central source may not be the same in Seyfert 
2 galaxies (despite the predictions of the unification models 
in their simplest form). Our results are consistent with the 
idea that the central source in ESO103-G35 is seen close to 
face on and at a larger inclination in the other objects. In 
the latter case, the contribution from the reflection compo- 
nent and the disk iron line emission will be supressed and 
it will be more difficult to detect any variations associated 
with them. 

The main result from the RXTE observations of the 
four Seyfert 2 galaxies ESO103-G35, IC5063, NGC4507 and 
NGC7172 is that they all show significant spectral variabil- 
ity in the sense that the spectrum steepens with increasing 
flux. This is a model independent result and is based on the 
analysis of the hardness ratios (see Section 4.1, Fig. 2). The 
origin of this variability may not be the same for all of them. 
In NGC4507 and NGC7172 we observe flux related N H and 
r variations respectively. In ESO103-G35 the spectral vari- 
ability could be due either to photon index variations (like 
NGC7172) or to reflection component changes. The results 
for IC5063 are puzzling. The stability of the iron line flux 



does not support reflection component variations as the rea- 
son of the observed spectral variations. On the other hand, 
the photon index variations are not well correlated with the 
continuum flux. Since this object is a radio galaxy, perhaps 
the normal "Seyfert - type" X-ray source in this case is "con- 
taminated" by emission from the non-thermal source which 
is responsible for the radio emission. 

The differences in the origin of the spectral variability 
observed in these objects suggest that the orientation of the 
nucleus and/or the distribution of the absorbing matterial 
around the nucleus is not the same in all Seyfert 2 galax- 
ies. We believe that future monitoring observations of these 
objects with XMM mission which has both a large effective 
area and good spectral resolution are necessary in order to 
clarify the issue of spectral variability in Seyfert 2 galax- 
ies and provide stringent constraints on the geometry of the 
nucleus and the location of the cirumnuclear matter. 
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Table 5. Power-law and Gaussian line fits for ESO103-G35 



Model 


N 5 


r 


Energy (keV) 


cr (keV) 


X 2 /dof 


A 1 
B 2 
C 3 
D 4 


14 l +0 -\ 
12.6+^-^.0+1° 
13.9l°' 4 
13.8+° j 


-, ,-7+0.02 

L - 01 -0.03 

1 rc + 0.03 

i - OO -0.03 
1 47+0.17 _ 1 61 +0.04 
± -^'-0.10 1 - ol -0.04 

1 K4+0.02 


5 QQ+0 01 

°' y,5 -0.01 

5 9i+° 04 

J - y -0.04 

5 90+ ' 04 

d - yu _0.04 
°- yu -0.05 


fiq+°'° 5 
u - Dy -0.04 

72+°-° 5 
U -' Z -0.05 
n 7 o+0.05 

0.1 


791.3/730 
684.3/716 
629.4/716 
723.7/716 



Notes: 1 N H , T tied; 2 N H untied, T tied; 3 Af H tied, T untied; 
4 gaussian line normalizations untied ; 5 in units of 10 22 cm~ 2 



Table 6. Power-law and Gaussian line fits for IC5063 



Model 




r 


Energy (keV) 


cr (keV) 


X 2 /dof 


A 1 
B 2 
C 3 
D 4 


19.2+1 7 
15.0+1° -28.0+ 2 ° 

n.9+_\% 

17.8l 2 j 


i c-y+0.07 

i cn + 0.10 

1-0 -0.09 
i ic+0.10 i 7n +0.20 
1 - 1O -0.13 1 -' U -0.19 

1 47+0.15 

1 -^'-0.06 


6 11+ - 08 

6 02+° 10 
D ' uz -0.03 

fi 05+°' 09 

6 04+° 13 
°- u ^-0.14 


07+0.13 

u - a '-0.17 

51+° 14 
U - Oi -0.06 

n 5n+°' 15 

u - ou _0.14 
5 o+0.09 


310.0/240 
191.9/236 
177.2/236 
302.6/236 



Notes: 1 N H , T tied; 2 N H untied, T tied; 3 N H tied, T untied; 
4 gaussian line normalizations untied ; 5 in units of 10 22 cm~ 2 



Table 7. Power-law and Gaussian line fits for NGC4507 



Model 


N 5 


r 


Energy (keV) 


cr (keV) 


X 2 /dof 


A 1 
B 2 
C 3 
D 4 


45.011J 
40.311° -50.9l 2 ;° 
45.2+2 
46.711 


1 40+0- 05 
1 -^ u -0.01 
1 40+0-05 

-, 97 +0.06 i co+0.05 
1 - z '-0.01 1 - oa -0.04 

1 49+0.03 

1 -^-0.03 


6 09+"- 04 
°' uy -0.06 

6 03+°' 05 
fi 05+°' 05 

5 qfi+0 03 
°- yD -0.03 


33+0.04 

u - oo -0.06 

42+ - 05 
U - 4Z -0.05 

40+°'° 4 
u -* u _0.03 

52+ ' 04 
u - oz -0.04 


887.2/807 
742.5/794 
763.7/794 
815.8/794 



Notes: 1 N H , T tied; 2 N H untied, T tied; 3 N H tied, T untied; 
4 gaussian line normalizations untied ; 5 in units of 10 22 cm -2 



Table 8. Power-law and Gaussian line fits for NGC7172 



Model 


N 5 


r 


Energy (keV) 


cr (keV) 


X 2 /dof 


A 1 
B 2 
C 3 
D 4 


n 5+°' 3 

ii -°-0.8 

11 o +1 ° - 14 +1 ° 

1J - -1.0 ly: ' u -1.0 

11 5+°' 7 
H.O_0.8 
ii O+0.5 
i4 - d -0.5 


i Qi+0.03 
1 - Si -0.05 
i 01+O.O6 
L - SL -0.05 
1 66+ 14 - 1 96+ ' 03 

1 - 8U -0.04 


6 19+ (M)9 

u - ±a -0.09 
fi 18+0.08 

°- 18 -o.io 
6 19+o.os 
u - ly -0.09 

D - ±( -0.08 


0.1 
0.1 
0.1 
0.1 


304.2/437 
274.6/429 
274.7/429 
299.8/429 



Notes: 1 N H , T tied; 2 N H untied, T tied; 3 N H tied, T untied; 
4 gaussian line normalizations untied ; 5 in units of 10 22 cm~ 2 



Table 9. Reflection component fits for ESO103-G35 



Model Njj T Energy (keV) cr (keV) X 2 /dof 



A' 2 


17 2+°' 3 


1.9 


6 07+°- 03 
D ' u '-0.03 
c nc+0.05 
6 - 06 -0.05 


n .+0.10 

u -°-o.io 


627.5/730 


B' 3 


15.911° -°19.011;° 


1.9 


n c+o.09 

u -°-0.03 


556.2/716 


C' 4 


17.5l°; 2 


1 R4+0- 08 9 n9+ - 05 
1 - S4 -0.08 ~ 2 - U2 -0.05 


6 07+ 06 

°- u '-0.06 


47+0.O8 

u -^'-0.04 


523.5/715 


D' 5 


17 9+O.2 
1( - z -0.2 


1.9 


6 06+°' 05 
D - UD -0.05 


51+ 08 
U-01_o.04 


604.7/716 


E' 6 


17 0+O.2 
1(.0_ Q 2 


1.9 


6 06+°' 04 
D - UD -0.05 


5i+o- 08 


525.8/716 



Notes: 1 in units of 10 22 cm~ 2 ; 2 N H , T tied; 3 N H untied, T tied; 4 N H tied, T untied; 
5 gaussian line normalizations untied ; 6 reflection component normalisation untied 
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Table 10. Reflection component fits for IC5063 



Model 




r 


Energy (keV) 


ex (keV) 


x 2 /dof 


A' 2 
B' 3 

c ,4 

D' 5 
E'6 


21.911° 
18.411° -31.01^ 

no 7+O.8 
ZZ *'-0.8 

21.911;° 
22.311;", 


1.9 
1.9 

-, ,,-+0.10 9 9n +0.15 
J" ' -0.10 _ ^- ZU -0.15 

1.9 
1.9 


fi ia+0-10 
b - 18 -o.io 
6 17+ - 05 

011 -0.05 
fi 1S+0-05 
D-1S-0.05 

6-201°;^ 

6-18±g;lg 


0.1 
0.1 
0.1 
0.1 
0.1 


269.3/241 
186.8/237 
181.8/236 
264.3/237 
176.7/237 



Notes: 1 in units of 10 22 cm~ 2 ; 2 N H , T tied; 3 N H untied, T tied; 4 N H tied, T untied; 
5 gaussian line normalizations untied ; 6 reflection component normalisation untied 



Table 11. Reflection component fits for NGC4507 



Model 


N 1 


r Energy (keV) 


cr (keV) 


X 2 /dof 


A' 2 
B' 3 
C' 4 
D' 5 
E'6 


50.4lJ;i 
47.9l°;«-53.4ll; 2 
51.911° 
50.811;° 
50.811;° 


1.9 6.14l°;°3 
1.9 6.13l°;° 4 

1 qn+°-° 5 _ 9 01+ 07 fi 13+0.04 
1 - au -0.05 z - ul -0.07 D - ±J -0.04 
1 9 6 13+ - 04 

1 9 6 14+0' 03 


0.1 
0.1 
0.1 
0.1 
0.1 


737.2/807 
701.5/795 
700.7/794 
732.0/795 
717.8/795 


Notes: 


1 in units of 10 cm - 


~ 2 ; 2 N H , T tied; 3 N H untied, T tied; 4 


Nh tied, r untied; 



5 gaussian line normalizations untied ; 6 reflection component normalisation untied 
Table 12. Reflection component fits for NGC7172 



Model 


n 1 r 

H 1 


Energy (keV) 


a (keV) 


X 2 /dof 


A' 2 
B' 3 

c ,4 

D' 5 
E'6 


12.2l°; 4 1.9 
9.9+1° - 14 - 0+ 3 o 19 

11 7+0.8 ' -, 74 +0.15 , qn+0.08 
1 - L ''-0.4 -0.22 - L - aJ -0.07 
12.2+0' 4 1.9 

12.ll°-j 1.9 


6 21+°' 09 

6 21+ 010 

6 20+°- 10 

6 20+°' 10 
D - zu _0.10 

6 20+°' 10 


0.1 
0.1 
0.1 
0.1 
0.1 


299.7/436 
273.5/428 
274.1/428 
297.1/428 
280.7/428 


Notes: 


1 in units of 10 22 cm" 2 ; 2 N H ,T tied; 3 N H 


untied, T tied; 


4 7V H tied, 


T untied; 



gaussian line normalizations untied ; 6 reflection component normalisation untied 
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